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Abstract 

Agricultural expansion and intensification are major threats to global biodiver- 
sity, ecological functions, and ecosystem services. The rapid expansion of oil 
palm in forested tropical landscapes is of particular concern given their high 
biodiversity. Identifying management approaches that maintain native species 
and associated ecological processes within oil palm plantations is therefore a 
priority. Riparian reserves are strips of forest retained alongside rivers in culti- 
vated areas, primarily for their positive hydrological impact. However, they can 
also support a range of forest-dependent species or ecosystem services. We sur- 
veyed communities of dung beetles and measured dung removal activity in an 
oil palm-dominated landscape in Sabah, Malaysian Borneo. The species rich- 
ness, diversity, and functional group richness of dung beetles in riparian 
reserves were significantly higher than in oil palm, but lower than in adjacent 
logged forests. The community composition of the riparian reserves was more 
similar to logged forest than oil palm. Despite the pronounced differences in 
biodiversity, we did not find significant differences in dung removal rates 
among land uses. We also found no evidence that riparian reserves enhance 
dung removal rates within surrounding oil palm. These results contrast previ- 
ous studies showing positive relationships between dung beetle species richness 
and dung removal in tropical forests. We found weak but significant positive 
relationships between riparian reserve width and dung beetle diversity, and 
between reserve vegetation complexity and dung beetle abundance, suggesting 
that these features may increase the conservation value of riparian reserves. Syn- 
thesis and applications: The similarity between riparian reserves and logged for- 
est demonstrates that retaining riparian reserves increases biodiversity within oil 
palm landscapes. However, the lack of correlation between dung beetle commu- 
nity characteristics and dung removal highlights the need for further research 
into spatial variation in biodiversity-ecosystem function relationships and how 
the results of such studies are affected by methodological choices. 



Introduction 

Agricultural expansion and intensification are currently 
among the main causes of decline in global biodiversity 
and ecosystem services (Phalan et al. 2013). However, 
large areas of agriculture will continue to be a key feature 
of our landscapes as the human population expands 
(Godfray et al. 2010). These cultivated landscapes can 
contribute to the persistence of biodiversity and delivery 



of ecosystem services, but appropriate, active management 
is required to achieve this (Garnett et al. 2013; Melo et al. 
2013). 

Successful management of biodiversity and ecological 
processes in tropical agricultural landscapes is especially 
important. Tropical landscapes are often particularly bio- 
diverse, highly productive for cultivation, and influence 
ecological functions and services on a global scale (Balm- 
ford and Whitten 2003). While primary forests are 
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critically important for conserving tropical biodiversity 
and ecosystem functions (Gibson et al. 2011) and 
once-logged forests in Southeast Asia also have high con- 
servation value (Edwards et al. 2011; Slade et al. 2011), 
focussing on these habitats alone is not sufficient. The 
area of land dedicated to crops or livestock is much 
greater than that in reserves or unmodified by humans 
(EUis and Ramankutty 2008), and the landscapes sur- 
rounding protected areas may strongly influence their 
success (Laurance et al. 2012). Moreover, modified agri- 
cultural landscapes can also be an important habitat in 
their own right (MendenhaU et al. 2012). 

Retaining areas of native vegetation along rivers can 
help maintain biodiversity and ecological functions within 
agricultural areas. These linear forest fragments are called 
riparian strips, buffer zones, stream management zones, 
or riparian reserves (the latter is used in Malaysia and 
hence in this paper). They are primarily retained because 
riparian forest reduces run-off into streams, improving 
water quality, and benefitting aquatic fauna (Sweeney 
et al. 2004; Mayer et al. 2007). Riparian reserves are also 
able to support forest-dependent communities of many 
terrestrial taxa, including birds, small mammals, and 
amphibians (Marczak et al. 2010). They are generally well 
protected legally and are a common feature of many agri- 
cultural landscapes (Lee et al. 2004; Barlow et al. 2010b), 
so offer a feasible, realistic option to improve biodiversity 
within cultivated areas. 

Nevertheless, the ecological roles of riparian reserves 
remain poorly understood, particularly in Southeast Asia. 
The majority of existing studies on riparian reserves focus 
on temperate regions, and particularly on bird species (see 
Marczak et al. 2010 for a review). To our knowledge, there 
are only 15 studies in tropical regions that evaluate the 
ecological characteristics of existing riparian reserves, all of 
which focus on the neotropics or north-western Australia 
(HiU 1995; Laurance and Laurance 1999; de Lima and Gas- 
con 1999; Graham and Blake 2001; Galindo-Gonzalez and 
Sosa 2003; Harvey et al. 2006; Medina et al. 2007; Gillies 
and St. Clair 2008, 2010; Lees and Peres 2008; Barlow et al. 
2010b; Norris and Michalski 2010; Rodriguez-Mendoza 
and Pineda 2010; Gilhes et al. 2011; Viegas et al. 2014). 

Here, we investigate the ecological impact of riparian 
reserves in the oil palm plantations of Sabah, Malaysian 
Borneo. Palm oil is now the world's primary vegetable oil, 
a major biofuel feedstock, and a component in many 
household products (Fitzherbert et al. 2008). Global 
annual production of palm oil more than doubled between 
1970 and 2010, with over 80% of the total now produced 
by Malaysia and Indonesia (FAO 2014). Production in 
Southeast Asia is still increasing and oil palm plantations 
are also likely to expand in west Africa and Amazonia 
(Butler and Laurance 2010; Foster et al. 2011). Establish- 
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ing successful conservation strategies in oil palm areas 
therefore has implications for landscapes across the world. 

We chose dung beetles (Coleoptera: Scarabaeidae: 
Scarabaeinae) as a focal group to assess the conservation 
value of the riparian reserves. Dung beetles are described 
as a "high performance indicator" for tropical regions: 
their community metrics vary with habitat disturbance or 
fragmentation, they show congruency with several other 
taxa, and are low-cost to survey (Gardner et al. 2008; 
Nichols and Gardner 2011). Dung beetles also provide 
important ecological functions such as dung removal and 
bioturbation (Nichols et al. 2008), which are of wider sig- 
nificance to entire ecosystems. 

Here, we assess whether riparian reserves support dung 
beetle communities and dung removal rates characteristic 
of larger areas of forest. We examine how the structure of 
riparian reserves could be managed to improve the extent 
to which they retain communities similar to those in 
logged forest. We also ask whether riparian reserves 
enhance the provision of dung removal services to the 
surrounding oil palm areas. 

Methods 
Study sites 

AH study sites were located within a 600 km^ area around 
and including the Stability of Altered Forest Ecosystems 
(SAFE) project site in Sabah, Malaysian Borneo (117. 5°N, 
4.6°E). The area is a mixture of twice-logged lowland dip- 
terocarp rainforest, acacia, and oil palm plantations, in 
which palms were planted between 1998 and 2011. Fur- 
ther details are given in Ewers et al. (2011). 

We selected 23 focal sites along river banks. Seven were 
located in logged forest (areas of continuous forest at least 
500 ha in size), seven in areas of continuous oil palm with 
no riparian reserve adjacent to the river, and eight in areas 
of oil palm with a riparian reserve adjacent to the river 
(Fig. SI shows these sites). One site was located in Maliau 
Basin primary forest reserve (70 km from the SAFE pro- 
ject) as a reference point, but there were no other primary 
forest sites near enough to allow spatial interspersion of 
replicate primary forest sites. As all the remaining primary 
forest in Sabah is already protected (Reynolds et al. 2011), 
evaluating the ecological characteristics of large areas of 
logged forest versus a network of smaller forest strips is 
more informative for future conservation action. All our 
sites were separated by at least 1.5 km, and the riparian 
reserve sites were at least 900 m (mean distance 3.3 km, 
standard deviation (SD) = 2.5 km) from the logged forest 
boundary. Given that dung beetle movement within a 48- 
h period is thought to be less than 500 m (Roslin 2000; 
Larsen and Forsyth 2005), it is therefore unlikely that 
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dung beetles were drawn to the riparian reserve traps from 
the logged forest areas. 

At each site, we set up a sampling grid of 12 points, 
consisting of four transects perpendicular to the river 
(Fig. S2). Transects were 100 m apart, with sampling 
points at 0 m, 50 m and 100 m from the high water line. 
The spacing of the grid conforms to standard methods of 
dung beetle sampling (Larsen and Forsyth 2005). Due to 
variation in width of the riparian reserves (mean 49 m, 
SD = 30 m, referring to forest width on one side of the 
river), where the riparian reserve was narrow, some points 
in these grids fell in the surrounding oil palm area. 

Data collection and analysis 

AH data were collected between the end of February and 
the beginning of July 2011. Seasonal changes in the low- 
land dipterocarp forests of Borneo are very limited 
(Walsh and Newbery 1999; Kumagai et al. 2005), and 
these months all fall in the slightly drier half of the year 
(Hamer et al. 2005). 

All analyses were carried out in R (R Core Team 2012) 
using the packages vegan (Oksanen et al. 2013), lme4 
(Bates et al. 2012) and nlme (Pinheiro et al. 2013). 

Dung beetle community and land use 

Dung beetles were collected using pitfall traps baited with 
25 g of human dung. Human dung attracts a wide variety 
of species (Davis et al. 2001; Larsen et al. 2006) and is 
recommended as a standardized bait in tropical forests 
(Marsh et al. 2013). Each trap consisted of a plastic cup 
(8-cm top diameter, 5.5-cm bottom diameter, and 12.5- 
cm depth) half-filled with a solution of water, detergent, 
and salt. The traps were protected from the rain with a 
cover and collected after 48 h. The order of sites was ran- 
domized, and traps were set at no more than two sites in 
each 48-h period. 

We could not obtain sufficient human dung to supply 
both the traps and the dung piles, so we used cattle dung 
for the dung removal experiment. Preliminary work in 
similar forest sites in Sabah shows that large cattle dung 
baits attract a similar species composition to smaller 
human dung baits, with the exception of some carrion 
feeding species found in higher abundances in human 
dung (Slade et al. 2011, E. Slade and D. Mann, unpubl. 
data). To compare species and dung removal results, we 
removed data on these carrion feeding species (n = 13, 
highlighted in Table SI) from all analyses apart from 
those testing for the effect of riparian reserve structural 
features on the entire dung beetle community. 

For each sampling point (trap), we calculated dung bee- 
tle abundance, the number of functional groups present 
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(using classifications based on diurnal vs. nocturnal activ- 
ity, body length, and method of dung removal after Slade 
et al. (2007)), a diversity (Shannon index), and total bio- 
mass. We weighed beetles from 24 species taken from 
across the whole range of body sizes (between 7 and 51 
individuals per species, average = 27, SD = 8) and used a 
polynomial regression to estimate biomass for the remain- 
ing species (Logio(mass) = —1.64 + 5.61*Logio(length) 
- 4.39*Logio(length)^ + 1.99*Loglo(length)^ = 0.982). 

For each site, we calculated /? diversity (mean Soren- 
sen's similarity index) and species richness (using cover- 
age-based rarefaction methods (Chao and Jost 2012) 
through the iNEXT online software (Hsieh et al. 2013). 
Coverage-based methods of rarefaction provide a more 
informative comparison of richness among multiple sam- 
ples than individual or sample-based methods of rarefac- 
tion as the ratio of species richness is not compressed 
(Chao and Jost 2012). Rarefied species richness could not 
be calculated at the trap level due to four traps having 
only one or two beetles. 

Wherever possible, we retained data at the highest spa- 
tial resolution (trap level) for analyses. For response vari- 
ables where this was the case (abundance, functional 
group richness, diversity, and biomass), we analyzed the 
effect of land use (logged forest, riparian reserve, or oil 
palm) with generalized linear mixed models, using tran- 
sect nested within site as a random factor. Where 
response variables could only be calculated at the site 
level iP diversity and rarefied species richness), we ana- 
lyzed the effect of land use with a generalized least 
squares model. For all models, appropriate error distribu- 
tions were specified and transformations or weight struc- 
tures (varldent function as described by Zuur et al. 
(2009)) applied where necessary. For all analyses testing 
for an effect of land use, we excluded data from points 
that fell outside of the forest strip at the narrowest ripar- 
ian reserve sites so that we were carrying out a true test 
for differences between the three land uses. Some traps 
were lost due to flooding or other disturbances, so data 
were only obtained from 201 traps in total (82 from 
logged forest sites, 43 from inside riparian reserves, and 
76 from oil palm sites). 

Differences in community composition across land uses 
were explored using de-trended correspondence analysis 
(DCA, vegan function "decorana"), which performs well 
as an ordination method for displaying similarity of trop- 
ical insect communities along an environmental gradient 
(Brehm and Fiedler 2004). We tested for significant dif- 
ferences in community composition using a permuta- 
tional analysis of variance (vegan function "adonis") with 
999 permutations and site as a grouping variable. 

To determine whether the relative abundance of the 
functional groups differed between logged forest and 
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riparian reserves, we ran a mixed model with abundance 
as a response variable and functional group, land use, and 
their interaction term as predictors. 

Dung removal and land use 

To record dung removal activity, uniform pats of 700 g 
of cow dung were set out at each sampling point (« = 12 
at each site) and collected after 24 h. Large herbi- 
vores, such as the tembadau or wild cow [Bos javanicus 
d'Alton), Asian elephant [Elephus maxitnus L.), and 
bearded pig (Sus barbatus MiiUer), occur within the study 
area so the experimental dung pats resemble those occur- 
ring naturally. Dung removal experiments were carried 
out at least 1 month after pitfalls traps were collected, in 
order to avoid interference but also remain close enough 
for dung beetles assemblages to be similar (Slade et al. 
2011). The order in which sites were visited was random- 
ized. The dung was frozen for a minimum of 24 h before 
the experiment to kill any invertebrates already present. 
Data on mass loss were corrected for evaporation using 
estimates from three evaporation controls set at each site. 
For the controls, the cow dung was placed in a flat-bot- 
tomed sieve with mosquito netting sealed around the top 
(both 1-mm mesh), to prevent entry of any dung beetles. 

The effect of land use on the mass of dung removed 
was analyzed with a general linear mixed model, with 
transect nested in site as random factors and a log trans- 
formation for the response variable. As with the data for 
beetle communities, for the riparian reserve sites, we only 
used data from within the forest strips (total n = 111: 84 
from logged forest sites, 44 from within riparian reserve 
vegetation, 84 from oil palm sites). 

We assessed whether the relationship between dung 
beetle community characteristics and dung removal was 
consistent across all land use types with a generalized 
linear model including land use, rarefied species richness 
(correlated with diversity = 0.6, P = 0.004), biomass 
(correlated with abundance, = 0.7, P = 0.0002), func- 
tional group richness, and all two-way interactions. As 
this analysis included coverage-based rarefied species rich- 
ness all other data were averaged to the site level. 

Dung beetle community structure and riparian 
reserve characteristics 

To analyse the effect of riparian reserve width and vegeta- 
tion complexity on the dung beetle community, we 
included data on all dung beetle species (both carrion 
and dung feeders). In order to test whether increasing the 
proportion of area in the riparian zone left as native 
vegetation impacts the dung beetle community, we 
included all the points in the sampling grid at each ripar- 
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ian reserve site {n = 95). Data were combined for each 
transect as this was the resolution at which riparian forest 
width could be measured (using GIS software [ArcMap 
version 10.1, ESRI, Redlands, CA]). A similar approach 
was used by Viegas et al. (2014) to test whether reserve 
width affects dung beetle communities in the Amazon. 

To assess the vegetation structure at each sampling 
point, we measured humus depth, canopy density (using 
a spherical densitometer), and basal area (using the angle 
point method (Bitterlich 1984)). We estimated the height 
of the tallest tree to the nearest 5 m using a ruler held at 
arm's length and a known reference height at the base of 
the tree. We scored the understorey vegetation density 
(below 2 m) and midstorey vegetation density (between 
2 m and 5 m) on an ordinal scale of sparse (fewer than 20 
stems or branches), medium (20-60 stems or branches), 
and dense (few patches of light and 60-100 + stems or 
branches). To obtain one numerical index summarizing 
the greatest variation in these data, we ran a metric scal- 
ing analysis on all these measurements. The first axis was 
positively correlated with canopy density, tree height, 
humus depth, basal area, and midstorey density. Because 
this output is therefore capturing variation in the three- 
dimensional structure of the habitats, we refer to it as a 
vegetation complexity index. 

We analyzed the effect of vegetation complexity on 
dung beetle abundance, biomass, diversity, functional 
group richness, and species richness using only data from 
sampling points falling within the riparian reserve forest. 
To test for any effects of reserve width or vegetation com- 
plexity, we used generalized linear mixed models, with 
site as a random factor and specified error families where 
appropriate. 

Provisioning of dung removal services by riparian 
reserves 

We analyzed the effect of riparian reserves on dung removal 
rates in the surrounding oil palm area in two ways. First, 
we compared the dung removed in oil palm adjacent to a 
riparian reserve (i.e., from sampling points at riparian 
reserve sites that fell outside the riparian forest, n = 52) 
and in oil palm without an adjacent riparian reserve (and 
also at least 50 m from the river bank, n = 56). Second, 
using only the data from sampling points in oil palm adja- 
cent to a riparian reserve, we analyzed the effect of distance 
fi-om the riparian reserve boundary on the mass of dung 
removed. We used a generalized linear mixed model with 
the presence/absence of riparian reserve or distance from 
the reserve boundary as a fixed factor for the two analyses, 
respectively. In both cases, we specified transect nested 
within site as a random factor and applied log transforma- 
tions to meet model assumptions. 
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Results 

Dung beetle community structure and land 
use 

In total, we identified 73 species from 9135 dung beetles 
(Table SI). The iNEXT software estimate for species 
coverage of the raw data was greater than 0.91 (>91% of 
species present were recorded) at all sites, and all spe- 
cies richness curves were approaching the asymptote 
(Fig. S3), indicating that we had sampled the community 
thoroughly. 



Table 1. Effects of land cover and habitat characteristics on dung 
beetle community metrics and dung removal using data at the trap 
level (or pooled to transect level for analyses with width as a fixed 
factor). Test statistics given for comparison of model specified against 
the null model (response - 1). 



Model 



df P 



Dung beetle community 

response to land use 
Abundance ~ land cover 
Biomass ~ land cover 
Shannon diversity ~ land cover 
Functional group count ~ land cover 

Dung beetle community 

response to reserve width 
Abundance ~ riparian reserve width 
Biomass ~ riparian reserve width 
Species richness ~ riparian reserve width 
Shannon diversity ~ riparian resen/e 

width 

Functional group richness ~ riparian 
resen/e width 

Dung beetle community 

response to vegetation complexity 
Abundance ~ vegetation complexity 
Biomass ~ vegetation complexity 
Species richness ~ vegetation complexity 
Shannon diversity ~ vegetation 

complexity 
Functional group richness ~ vegetation 

complexity 

Dung removal function 
Dung removed ~ land cover 
Dung removed in oil palm ~ the 

presence/absence of riparian reserve 
Dung removed in oil palm ~ distance 

to riparian reserve boundary 

Riparian reserve features 
Vegetation complexity ~ riparian 
resen/e width 



5.9 
7.9 
22.7 
28.8 



0.62 
0.05 
3.69 
5.45 



5.95 
0.54 
0.3 

0.0004 



0.0005 1 



4.6 
0.58 



2.11 



0.051 
0.019* 
<0.0001» 
<0.0001' 



0.43 
0.82 
0.055 
<0.02* 

0.28 



0.015* 
0.46 
0.58 
0.98 

0.98 



0.10 
0.45 



0.15 



0.1? 



Significant differences between the model described and the null 
model (*P < 0.05; **P < 0.01; ***P < 0.001). 



We found a significant effect of land use on dung bee- 
tle biomass and a weakly significant effect of land use on 
abundance (Table 1). The biomass in riparian reserves 
was intermediate between oil palm and logged forest, but 
not significantly different from either (Fig. IB). 

Species richness (Table 2), diversity, and functional 
group richness (Table 1) also varied significantly with 
land use. Riparian reserve species richness (Fig. IC), 
diversity (Fig. ID), and functional group richness 
(Fig. IE) were significantly lower than in logged forest 
and higher than in oil palm. We found no difference in 
within-site P diversity (mean Sorensen's similarity index) 
among the different land use types (Table 2). 

Three of the seven functional groups were missing 
completely from oil palm sites: small nocturnal tunellers, 
small diurnal rollers, and large nocturnal rollers. In con- 
trast, aU functional groups were found in at least one of 
the riparian reserves sites. However, we found a signifi- 
cant interaction between land cover and functional 
group on dung beetle abundance (x^ = 59.8, df = 6, 
P < 0.0001). This indicates that functional groups vary in 
the extent to which they are negatively impacted by the 
conversion from logged forest to riparian reserve 
(Table 3). The most negatively impacted functional 
groups were the large diurnal tunellers and large diurnal 
rollers. 

The community composition of the riparian reserves 
was more similar to logged forest than oQ palm, although 
a distinct difference in the communities of the reserves 
and larger forested areas remains (f 1,193 = 21.4, 
P = 0.001, Fig. 2). 

The single primary forest reference site (not included in 
the analyses above) had much higher mean dung beetle 
abundance (145% of logged forest), biomass (319% of 
logged forest), diversity (115% logged forest), and func- 
tional group richness (114% logged forest) than all other 
land use types (Fig. 1). However, the species richness of the 
reference primary forest site fell within the range of the 
logged forest sites (Fig. IC). 

Dung removal and land use 

The proportion of dung removed across all sites was low 
(mean = 0.1, SD = 0.14). There was no significant rela- 
tionship between dung removal and land cover, species 
richness, diversity, functional group richness, or any of 
the two-way interactions (Table 2). 

Dung beetle community and riparian 
reserve characteristics 

There was no significant relationship between riparian 
reserve width and vegetation complexity (Table 1). 
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Forest Riparian reserve Oil palm 



(F) o 

m 




Forest Riparian reserve Oil palm 



Figure 1. Effect of land cover on (A) dung 
beetle abundance, (B) biomass, (C) coverage- 
based rarefied species richness, (D) diversity 
(Shannon index) (E) functional group richness 
and (F) dung removal. All panels show means 
and standard errors. The dotted lines indicate 
values for the one primary forest reference site 
(for visual comparison only; the data were not 
included in the analysis). Stars denote 
significant differences between groups based 
on model contrasts (*P < 0.05, **P < 0.01, 
***P < 0.001). 



Table 2. Effects of land cover and habitat characteristics on dung 
beetle community metrics and dung removal on response variables 
that could only be calculated at site level. 



Model 



df 



Species richness ~ land cover 

fS diversity ~ land cover 

Dung removed ~ sp.rich*land.cov + 

biomass*land.cov + f.rich*land.cov 
Dung removed ~ sp.rich + biomass - 

f.rich 



16.9 2,19 <0.0001*** 
1.9 2,19 0.18 
0.8 11,10 0.67 



1.7 



3,18 0.2 



sp.rich, species richness; land.cov, land cover; f.rich, functional group 
richness. 

Significant differences between the model described and the null 
model (*P < O.OB; **P < 0.01; ***P < 0.001). 



We found no evidence of a relationship between 
reserve width and dung beetle abundance, biomass, 
or functional group richness (Table 1). However, there 
was a significant positive relationship between width 
and diversity and a weakly significant positive relation- 
ship between riparian reserve width and dung beetle 
rarefied species richness (Shannon index, Table 1, 
Fig. 3). 

We found a positive relationship between the vegeta- 
tion complexity of the riparian reserve forest and bee- 
tle abundance (Table 1). However, we found no 
significant effects of vegetation complexity on biomass, 
species richness, diversity, or functional group richness 
(Table 1). 
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Table 3. Model output of generalized linear mixed model (dung bee- 
tle abundance - land cover * functional group), showing parameter 
estimates and standard error for the percentage decline in abundance 
of each functional group in riparian reserve sites relative to logged 
forest sites. 





No. 


Estimate 


Standard 






species 


% 


error % 




Functional group 


in group 


decline 


decline 


P 


Large diurnal tunellers 


1 


89.48 


17.70 


0.0002*** 


Large diurnal rollers 


2 


80.19 


15.78 


0.0002*** 


Small diurnal rollers 


1 


71.92 


12.74 


0.0001*** 


Small diurnal tunellers 


21 


18.65 


5.60 


0.0054** 


Large nocturnal 


4 


32.96 


18.64 


0.3 


tunellers 










Small nocturnal 


3 


95.32 


71.8 


0.2 


tunellers 










Large nocturnal rollers 


1 


86.48 


207.1 


0.7 



Significant denotes (*P < 0.05; **P< 0.01; ***p< 0.001). 
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Figure 2. Detrended correspondence analysis (DCA) plot indicating 
that riparian reserve community composition is more similar to forest 
than oil palm. Ellipses show standard deviation around the mean for 
each land use. 

Provisioning of dung removal services by 
riparian reserves 

Dung removal did not differ significantly between oil 
palm with and without riparian reserves, and we found 
no significant effect of distance from the riparian reserve 
boundary on the mass of dung removed (Table 1). 

Discussion 

The rapid expansion of oil palm plantations throughout 
the tropics threatens many forest species. While large pro- 
tected areas will undoubtedly remain the priority for con- 
servation in these areas, riparian reserves are a potential 
opportunity to increase the biodiversity retained within 
agricultural areas. Our data show that, compared with 
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areas of oil palm, riparian reserves of at least 30 m width 
(on each side of the river) support dung beetle communi- 
ties more similar to those in adjacent areas of logged 
forest. We found some effects of reserve width and vege- 
tation complexity on the diversity and abundance of dung 
beetles in riparian vegetation, suggesting that these struc- 
tural features may also make a limited contribution the 
biodiversity benefits of the riparian reserves. 

Previous studies on dung beetle communities and dung 
removal function in Bornean rainforests have shown a 
positive relationship between dung beetle species or func- 
tional group richness and dung removal function (Slade 
et al. 2007, 2011). Our data suggest that this relationship 
may not hold in riparian zones and that the presence of a 
riparian reserve does not increase dung removal in sur- 
rounding areas of oil palm. 

Dung beetle community structure and land 
use 

We found that dung beetle communities in oil palm plan- 
tations had lower biomass, species richness, diversity, and 
functional group richness than larger areas of forest. Sim- 
ilar studies in Borneo have also found that dung beetle 
species richness and diversity declines with logging and 
conversion to plantations (Davis et al. 2001; Edwards 
et al. 2013). In contrast to these studies, we did not find 
a decline in dung beetle abundance across the land use 
gradient, but this may be because we did not make a 
comparison with primary forest sites. 

The species richness, functional group richness, diver- 
sity, and overall community composition of dung beetle 
communities within riparian reserves were more similar 
to forest than oil palm. These results provide strong evi- 
dence that protecting riparian reserves retains biodiversity 
within oil palm landscapes, when compared to plantations 
where oil palm is planted up to the river bank. Riparian 
reserves are not, however, an adequate substitute for large 
areas of logged or primary forest. In addition, farther 
work is needed to establish whether dung beetles are able 
to maintain viable, self-supporting populations within 
riparian reserves, rather than comprising transient visitors 
or sink populations (Barlow et al. 2010a). 

The functional diversity of the oil palm was lower than 
both the logged forest and riparian reserves. Rollers and 
nocturnal species were particularly negatively affected by 
conversion to oil palm. The greater sensitivity of these 
groups after conversion to oil palm was also reported by 
Edwards et al. (2013) and may be connected to limited 
temperature tolerance. The riparian reserve network 
retained all functional groups, but each site tended to 
support fewer functional groups than logged forest sites. 
The drop in average functional group richness between 
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Figure 3. Relationship between riparian reserve width and (A) rarefied species riclnness and (B) diversity (Sliannon index). Plots sinow 
mean ± standard error for eacin replicate site. 



forest and riparian reserves was due to the loss of diurnal 
species (both small and large, rollers and tunellers), which 
may be due to a decline in diurnal mammal species 
(Andresen and Laurance 2007). These results suggest that 
the isolation of forest strips results in different trait- 
dependent responses compared to conversion to oil palm. 
Our findings contrast the global study carried out by Nic- 
hols et al. (2013a), who found that for the afro-eurasian 
tropics nocturnal species are more affected by forest mod- 
ification, and diurnal species are more negatively affected 
by conversion to plantation. However, Nichols et al. 
(2013a) were comparing all forest modification to a pri- 
mary forest baseline, whereas we are comparing oil palm 
and riparian reserves to a logged forest baseline, which 
may explain this discrepancy. 

As well as supporting dung beetle species that would 
not survive if oil palm were planted along river banks, 
riparian reserves are likely to benefit a range of other 
taxa. In Borneo, riverine forest corridors are recognized 
as important habitat for some mammalian species, 
including the orang-utan, proboscis monkey, and pygmy 
elephant (Venkataraman et al. 2009), but little research 
has been carried out on the importance of riparian 
reserves for many other groups in this region. Riparian 
reserves in the neotropics support communities of birds, 
amphibians, and small mammals found in undisturbed 
forest (de Lima and Gascon 1999; Lees and Peres 2008) 
and also facilitate movement of forest specialists through 
agricultural land (Gillies and St. Clair 2008). Alongside 
the hydrological benefit of riparian reserves, their role 
in conserving terrestrial species should be more widely 



recognized by sustainable management guidelines. This is 
especially the case in Sabah where all the remaining pri- 
mary forest is already protected and increasing conserva- 
tion in cultivated landscapes is arguably the highest 
priority (Reynolds et al. 2011). 

Dung removal and land use 

Despite the significant differences among beetle commu- 
nities in different land uses, we did not detect a signifi- 
cant effect of land cover on dung removal rates over 
24 h, nor any significant relationship between dung 
removal and species richness, biomass or functional group 
richness. This contrasts a number of studies that show 
strong positive correlations between dung beetle species 
richness or functional group richness and dung removal 
rates (Slade et al. 2011; Braga et al. 2013; GoUan et al. 
2013), and also with evidence that dung removal rates in 
Amazonian riparian reserves are higher than in surround- 
ing pasture (Norris and Michalski 2010). 

There are several possible explanations for these results. 
First, dung removal rates were low (e.g., compared with 
those in primary forest nearby; Fig. IF), and it is possible 
that a difference in removal would be seen if dung pats 
were left out for longer, and a greater proportion of mass 
was removed. Secondly, it is possible that differences in 
the communities attending the two bait types diminish 
our ability to detect correlations between biodiversity and 
function; Nichols et al. (2013b) discuss how dissimilarities 
between the response of dung beetles communities and 
dung burial rates to human impact in the Amazon may 
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be an artifact of surveying the community and function at 
separate times with different baits. However, correlations 
have previously been detected using different baits within 
similar forests in Borneo (Slade et al. 2011). Thirdly, it is 
possible that there is spatial variation in the relationship 
between dung beetle community composition and dung 
removal. Because the mortality of dung beetle larvae may 
increase with soil moisture content (Sowig 1995), dung 
beetles may not build nests (bury dung) near rivers even 
though they come to baits to feed. The positive relation- 
ship between species richness and dung removal may 
therefore break down in riparian zones, but additional 
data on how dung removal and soil moisture vary is 
needed to confirm this. Therefore, while there may be 
local and regional variation in biodiversity-ecosystem 
function relationships within tropical forests, the extent 
to which these relationships are affected by sampling 
methodology needs to be further resolved. 

The lower functional group richness in riparian reserves 
compared with forest may affect important ecological 
processes inside the reserves. In particular, roller species 
are less abundant in the reserves. As these species roll 
dung balls, often containing seeds, horizontally away from 
the dung pat they can potentially reduce the negative 
effects of seed clumping and seedling competition (Law- 
son et al. 2012). It is therefore possible that germination 
and dispersal dynamics of plant species in the reserves are 
impaired relative to logged forest. Other processes such as 
soil bioturbation, soil fertilization, and parasite suppres- 
sion that are mediated by dung beetles (Nichols et al. 
2008) may also be reduced in the riparian reserves as a 
result of the decline in some functional groups. 

Dung beetle community structure and 
riparian reserve characteristics 

Our results suggest that the width and vegetation complex- 
ity of riparian reserves may have a positive impact on dung 
beetle diversity and abundance, respectively. Although 
these relationships are weak, these findings are of direct 
relevance to management and policy specifications. Legal 
requirements for the protection of riparian forest exist in 
a number of countries (Sabah Water Resources Enact- 
ment, 1998; Barlow et al. 2010b; Marczak et al. 2010; 
FAO, 2014), and riparian reserves are also included in the 
criteria for certification of sustainable palm oil produc- 
tion (e.g., by the Round Table on Sustainable Palm Oil 
(Criterion 4.4 (RSPO 2013)). However, very little ecologi- 
cal information has influenced the details of these guide- 
lines (Barlow et al. 2010b; Ewers et al 2011). Our 
findings highlight the need for further research to clarify 
the importance of the structural features of riparian 
reserves. 
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Provisioning of dung removal services by 
riparian reserves 

Grazing of cattle underneath oil palms is expanding in 
Malaysia (Latif and Mamat 2002), and the requirement 
for dung removal services within these landscapes is likely 
to increase. However, our results suggest that retaining 
riparian reserves within oil palm plantations may not 
contribute to an increase in dung removal services within 
surrounding oil palm. 

Conclusions 

Overall, it is evident that riparian reserves can contribute 
toward the conservation of dung beetle communities that 
are threatened by the expansion of oil palm, but that the 
extent to which they support dung removal activity and 
other terrestrial ecosystem services requires further study 
over greater spatial and temporal scales. Nevertheless, the 
results presented here indicate that riparian reserves 
should be more widely recognized as a conservation strat- 
egy for terrestrial biodiversity. We must emphasize that 
we do not recommend riparian reserves within oil palm 
plantations as an alternative to protecting large areas of 
primary or secondary forest. On the contrary, we feel that 
as an addition to such protection, riparian reserves should 
be more widely recognized as a promising opportunity 
for conservation in tropical agricultural landscapes. 
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